Band Structure, Phase transitions and Semiconductor Analogs in One-Dimensional 

Solid Light Systems 
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The conjunction of atom-cavity physics and photonic structures ("solid light" systems) offers 
new opportunities in terms of more device functionality and the probing of designed emulators of 
condensed matter systems. By analogy to the canonical one-electron approximation of solid state 
physics, we propose a one-polariton approximation to study these systems. Using this approximation 
we apply Bloch states to the uniformly tuned Jaynes-Cummings-Hubbard model to analytically 
determine the energy band structure. By analyzing the response of the band structure to local 
atom-cavity control we explore its application as a quantum simulator and show phase transition 
features absent in mean field theory. Using this novel approach for solid light systems we extend the 
analysis to include detuning impurities to show the solid light analogy of the semiconductor. This 
investigation also shows new features with no semiconductor analog. 

PACS numbers: 42.50.Pq, 64.70.Tg, 71.36.+C 



I. INTRODUCTION 



Quantum properties of strongly correlated many- 
particle systems are intensive to compute and experimen- 
tally hard to access in conventional solid state materi- 
als. Because of this, artificial tunable systems to sim- 
ulate highly correlated many-body dynamics have been 
proposed. Early attempts at these so called quantum 
simulators used arrays of Josephson junctions [l|, 0] to 
reproduce strong many-body bosonic interactions. Cold 
atoms trapped in optical lattices Q now constitute an im- 
portant experimental platform for quantum simulation. 
Recently, arrays of coupled atom-photon cavities have 
been proposed 0, E E] as another highly flexible plat- 
form for quantum simulation. This has opened up excit- 
ing new possibilities for combining atom-cavity and pho- 
tonic structures. A two-level atom in an optical cavity is 
described by the well known Jaynes-Cummings (JC) Q 
model that gives rise to nonlinear photon- pho ton inter- 
actions including photon blockade [1, E [Tolllll| . By con- 
sidering a network of evanescently-coupled JC systems, 
a tight binding model can be developed, which is the 
Jaynes-Cummings-Hubbard (JCH) model [IEEE- 
In this paper, we consider this particular model and by 
determination of its band structure in response to local 
atom-cavity control we explore the application of the sys- 
tem as a quantum simulator. 

A number of realizations of the JCH model are possible 
including photonic bandgap (PBG) structures with sin- 
gle atoms Q or small ensembles [l4[ , arrays of supercon- 
ducting strip-line cavities II 91 . cold atoms connected via 
optical hbre interconnects [161 ] , plasmonics [17| , or linear 
ion traps [H, E- For clarity we focus on the concep- 




> ' ' V> ' ' ' ??? ' 



7 \H 



'Electronic address: jamesq@unimelb.edu.au 



FIG. 1: (Color online) A possible realisation of a one di- 
mensional Jaynes-Cummings-Hubbard model. Here holes are 
drilled into a thin membrane and lattice defects serve as the 
optical cavities housing two-level atoms. The cavities are 
coupled evanescently with single-photon hopping rate k and 
the atom is coupled to the local cavity mode with coupling 
strength j3. 



tually simple PBG implementation. A photonic crystal 
is realized by a periodic modulation of refractive index 
that is arranged so as to have an optical bandgap i.e. a 
range of photon energies that can not propagate through 
the medium [2fj| |. Typical implementations are arrays 
of holes drilled in a high refractive index membrane [2 11 ] . 
Lattice defects, i.e. where holes are not drilled, act as op- 
tical cavities into which a two-level atom can be placed as 
illustrated in Fig. [T] If the membrane was fabricated in 
diamond, the two-state systems could be ion-implanted 
NV centers, and there exists suitable designs for such 
cavities [22, 23] . The photon-atom coupling provides the 
mechanism by which the polaritons strongly interact. 

An infinite JCH lattice exhibits a quantum phase tran- 
sition analogous to the insulator-superfluid transition of 
the Bose-Hubbard model. This has been confirmed by 
several methods including mean-field (MF) [g, |2J] , den- 
sity matrix renormalization group [251] . v ariational clus- 
ter [26| . and linked-cluster expansion [l3|, and quantum 
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Monte Carlo [27[. Zhao et. al [28[ have also argued that 
the JCH phase transition belongs to a different universal- 
ity class from the Bose- Hubbard model. The signature of 
transition has also been identified in finite cavity simula- 
tions |12|. The term "solid light" was so coined because 
of the presence of the Mott insulation (MI) phase in the 
systems. 

The macroscopic to mesoscopic size of conventional 
solid state systems, has meant that it has been natural 
in solid state physics to assume infinite lattices. Tra- 
ditional cavity QED systems have been restricted to few 
(usually one) cavities , ho wever this is changing due to ad- 
vances in fabrication [29( . These new opportunites reveal 
a parameter space where condensed-matter concepts can 
be usefully applied to quantum optical systems, which 
in return act as quantum simulators to study non-trivial 
condensed matter systems. 

Within this context, we consider the case of a one- 
dimensional (ID) infinite chain of cavities. By analogy 
with the canonical one-electron approximation in solid 
state physics, we discuss the one-polariton approximation 
in a solid light context. This work informs and extends 
the literature on ID coupled-cavity waveguide photon 
propagation [3(1 HH, H2] • In particular, following a review 
of the JCH model in SecHI] we apply Bloch states to an- 
alytically develop the band structure of solid light, and 
we ascertain the phase transition diagram of Mott lobes 
and the superfluid (SF) regions in Sec. IIIII One should 
note that, while the band structure of PBG systems is 
well understood, we are exploring the band structure of 
the defects within the PBG, i.e. bands within the band 
gap. In Sec. IIVI we extend the one-polariton analysis to 
investigate how semiconductor properties can be simu- 
lated in the JCH model by embedding detuned atoms 
in a tuned system, showing a solid light analog of the 
extrinsic semiconductor. 



II. THE JAYNES-CUMMINGS-HUBBARD 
MODEL 

The Jaynes-Cummings (JC) model describes a two- 
level atom interacting with a quantized mode of an opti- 
cal cavity. It is described by the Hamiltonian (H = 1), 

H JC = ea+o-- +wo t a + (3(d + a + ^a t ), (1) 

where cr + , o~ (a\ a) are the atomic(photonic) raising and 
lowering operators respectively, e is the atomic transition 
energy, u> is the cavity resonance frequency, and (3 is the 
atom-photon coupling. 

The eigenvectors of the JC Hamiltonian are the dressed 
states (polaritons), 

|±,n) = f=== ,V n > 0, 



with corresponding eigenvalues, 

E\±,n) =nuj± x(n) - A/2, 



(2) 
(3) 



where \g,n) and \e,n — 1) are the atomic ground and 
excited states with n excitations, detuning parameter 
A = u> — e, and n photon generalized Rabi frequency 
x(n) = \/n/3 2 + A 2 /4. For n — 0, the eigenstate is \g, 0) 
with zero eigenvalue. 

A feature of the JC model is photonic blockade, 
whereby absorption of a photon blocks absorption of a 
subsequent photon because the required transition is de- 
tuned from resonance. This effect provides a mechanism 
by which polaritons can strongly interact. Strong atom- 
photon coupling is required for photonic-blockade to be 
observed and ideally very small cavities with very high 
quality factors, Q, are necessary. 

To go beyond one-cavity systems, we consider the Hub- 
bard model [lH which is one of the simplest of the models 
describing a system of interacting particles in a lattice. 
The Hubbard model parameterizes the Hamiltonian by 
two energy scales, those of inter-site hopping, and on- 
site repulsion. Quantum phases are seen as the rela- 
tive strength of these interactions are varied at a defined 
chemical potential. 

Combining the JC and Hubbard models, the JCH 
model describes an array of Jaynes-Cummings (JC) cavi- 
ties with Hubbard-like inter-cavity interaction where the 
JC interactions provides an effective particle-particle in- 
teraction. Grouping the JCH Hamiltonian into intra and 
inter-cavity terms yields, 



H 



E^ JC 



(r,s) 



.ala. 



E 



H r N r 



(4) 



where r and s are cavity site indices, [i is the chemical 
potential, and N is the total number of cavity excitations. 

The chemical potential energy is the change in mean 
energy per additional particle, /x = dE/dN . While this 
is typically zero in photonic systems, the nonlinear inter- 
actio ns g enerated by the JC Hamiltonian makes fi non- 
zero [l2l l3~ij . One should realize that in our case the 
chemical potential is not an independent parameter. In 
general the cost of adding particles is a function of all of 
the system parameters, in particular the strength of the 
atom-photon interactions and the number of excitations. 
Experimentally, /i can be determined via spectroscopy 
analysis of the change in system energy with excitations. 
However, for easy comparison with the standard phase 
transition literature, we here treat \i as an independent 
parameter. 



III. THE ONE POLARITON APPROXIMATION 
AND BAND STRUCTURE 

Because of the strongly interacting nature of the po- 
laritons, the JCH system becomes a complex many-body 
problem, akin to that of the electronic many-body prob- 
lem of solid state physics. Complexity quickly increases 
in difficulty with number of particles making exact so- 
lutions difficult to obtain. An approximation which has 
proven fruitful in describing the electronic properties of 
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solids is the so called one-electron approximation [35[ . In 
this formulation one decouples a single electron from the 
rest, treating the others as a mean sea of electrons. Thus 
along with the potential provided by the fixed lattice the 
single electron also sees a mean potential provided by the 
sea of charge distribution of the other electrons. Conse- 
quently, the total state of the system can be expressed as 
a linear combination of the states of one electron. 

Motivated by the one-electron approximation we inves- 
tigate the JCH system in a one-polariton approximation, 
describing an extra polariton above a uniformly filled lat- 
tice of JC systems. In the photon-blockade regime, we 
can also consider the concept of a polariton hole, i.e. a 
cavity with one fewer polariton than in the rest of the 
lattice. The one-polariton (hole) approximation is most 
valid in the Mott phase and, as we will show, captures 
the transition to the superfluid phase, although the ap- 
proximation will not be valid to describe the SF phase. 
The one-polariton (hole) approximation is defined as a 
reduced basis set, 



(5) 



where r and s indicate atomic sites, (p l m represents a site 
state with m total atomic and photonic excitations, I 
is indexed over all unique states of m excitations. The 
repeated index I is summed over all states with m ex- 
citations. For example in the bare basis, the set of I 
would be the ground and excited atomic states. The 
state co-efficients are such that the normalization condi- 
tion J^i c f — 1 is satisfied. For the basis set to be most 
valid in the Mott phase we choose q such that in the limit 
of zero hopping potential (i.e. k = 0), each site is a JC 
eigenstate. In the bare basis c; are just the co-efficients 
in the dressed states, i.e. 



c g (n) 
c e (n) 



^2x 2 (n)Tx(n)A 

-A/2± X (n) 
y / 2 X 2 (n) TX(«)A 



(6) 
(7) 



For notational convenience, in the rest of this paper 
we will adopt the one-electron convention where the s 
identical sites of the basis set are implied. 

We restrict our attention to the low-temperature limit 
of a ID infinite periodic chain of uniform sites, as de- 
picted in Fig. [Tj The onsite energy is given by, 



^H rs \(f>,s) = E r 



(8) 



where H rs is the Hamiltonian that relates site r to site s 
and E r is the onsite energy. 

The one-polariton approximation allows us to use the 
usual ansatz of Bloch's theorem for periodic structures 
to solve for the onsite eigenenergy. The ansatz is the 
discretized Bloch state, 



where d s is the displacement to site s, and k is the 
wavevector associated with the crystal momentum. We 
normalize units so that distance between neighboring 
sites is unity and k is in units of \d s — dr\ . 
Applying the Bloch ansatz to Eq. © we get, 



^i7 rs exp[ifc • (d s - d r )]\4>,r) = E r \</>,r) 



(10) 



which for the ID lattice becomes, 

[fr r , r + 2fl r , r +icos(A!)]|^r) = £7 P |0 1 r) , (11) 



where H r r = — I^N r is the onsite Hamiltonian and 

H r .r+i = — k(gv°v+i + a r+i c 4) is the hopping Hamilto- 
nian. As we have assumed identical sites, the chemical 
potential ji and hopping potential k are site independent 
and constant. 

The (2n x 2n) matrix of Eq. (Ilip is block diagonal, 
and each (2 x 2) diagonal block is the same function of 
n. The eigenenergy of the extra polariton is (n > 0), 



E±(n, k) = lo — fx + \(n) — h(n, k) 



± \A 2 (™,fc)-£K™,fc) + x(" + i) 2 , 



(12) 



where, 



h(n,k)= [n + (? g (n)}n cos(fc) , (13) 
g(n, k) = {A(n + l)y/n(3c e (n)c g (n) 

+ [(n + l)c g (n) 2 - nc e (n) 2 ] A}/tcos(fc) . (14) 

The energy of the hole is (n > 0), 

E^(n, k) = + x(n) - h'(n, k) 



± \fh' 2 (n, k) - g'(n, k) + x(n - If 



(15) 



where, 



, r) exp[ifc • (d s - d r )], 



(9) 



h'(n,k)= [n- c 2 e (n)]KCOs(k) , (16) 
g'(n, k) = {4(n - l)^(3c e (n)c g (n) (17) 
+ [(nc g (n) 2 — (n — l)c e (n) 2 ] A}k cos(fc) . 

Fig.[2]shows the one-polariton(hole) band structure. It 
forms the allowed energies a polariton can take near the 
Mott phase, analogous with conventional band structures 
which give the allowed energies an electron can take. 

An energy gap at zero temperature indicates the pres- 
ence of a MI phase. When the energy gap is zero, it 
costs no energy to add a polariton or hole, and we are 
no longer in the MI phase. This transition border oc- 
curs at critical chemical potential (i%(n, 0) = E v _ (n, 0) +[i 
and jii^(n, 0) = £^(n,0) - fj,. fi^(n, 0) will form the up- 
per boundary of the nth Mott lobe, and /xjf(n, 0) will 
form its lower boundary. The exception is n — 0, where 
there is no lower boundary. Fig. [3Ja)-(c) shows the band 
structure and critical chemical potential spectrum at var- 
ious points, and phase transition diagram. Overlayed in 



FIG. 2: The solid (dashed) lines are the one polariton (hole) 
band structures for A = 0, (ji-w)//3 = -0.5, n/f) = 0.01. 
The numbers indicate the n value. 



(b) 





FIG. 3: (Color online) (a) The dashed (solid) lines plot the 
one polariton- hole band structure at point A (A'). At A (A') 
the energy cost of adding a hole is less (more) than the po- 
lariton counterpart, (b) The dashed (solid) lines plot fi^' h at 
point B (B'). The gap at k = gives the width of the Mott 
lobes, (c) Plot of the energy gap indicating the Mott lobes. 
Overlayed is the phase transition boundary. Indicated also 
are the Mott lobe states. 



Fig. El^c) is the aforementioned phase transition bound- 
ary. 

Fig. GHa) gives the polariton and hole energy bands 
for point A and A' in Fig. [3jc). Just above the curve 
corresponding to /x^ (e.g. point A) the energy cost of 
adding a hole is less than that of adding a polariton. 
Just below the curve corresponding to yP c (e.g. point A') 



Kip 



FIG. 4: (Color online) Comparison of one-polariton approxi- 
mation (solid line), MF (dotted line), and finite cavity simu- 
lations (colored dashed lines). One mean excitation plateaux 
boundary are calculated for 2 to 5 cavities. Two mean exci- 
tiation plateaux boundary for 2 to 4 cavities. As the number 
of cavities increases it appears that finite cavity simulations 
tend to one-polariton analysis and not MF. 



the cost of adding a polariton is less than that of a hole. 
These minimal costs are the energy band gaps. Fig. [3jc) 
plots the energy gap forming the Mott lobes. 

The size of the gap indicates the level of stability, with 
larger gaps being more stable. The widths of the Mott 
lobes are given by (J%(n, 0) — n^(n,0). Fig. (3jb) plots 
(n, k) spectrum near and far from the lobe tip. 

As a means of comparison we overlay the one- 
polarition (hole) approximation phase transition border 
with that determined by a MF analysis in Fig. [4] In 
the MF analysis the decoupling approximation a\a r — 
(a! r )a r +{a s )a\— (a s ) (aj.) is made, so that the JCH Hamil- 
tonian is expressed as a sum over single sites, 



E 



[H 3 S C - ZKip(a\ + a s ) + zk\^\ 2 - fj,N s 



(18) 

where ip = (a) is known as the SF order parameter. It has 
been proposed as a measure of superfluidity [3||, where 
tp = indicates MI, and ip > indicates superfluidity, z 
is the number of neighbors. 

In Ref. [ll] it was shown that the lobes in the MF 
approximation correlate well with the plateaux of mean 
excitations seen in finite cavity systems, demonstrating 
the progression from the finite case to the infinite. Fig. [4] 
shows the connection between these finite calculations 
compared with the MF and one-polariton approximation 
about the first and second Mott lobe. As the number 
of cavities increases we see that the exact results for the 
quantum phase transition points tend towards the one- 
polariton approximation and not the MF approximation. 
This highlights the physics missed in the MF analysis, 
but captured in our treatment, namely the spatial super- 
position states of the single polariton. 
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IV. SEMICONDUCTOR SIMULATOR 

The solid light system with tunable onsite and inter- 
site coupling parameters and extensions to 2D and 3D 
is able in principle to capture the physics of many non- 
trivial condensed matter systems. As an introductory 
example we consider a case of the solid light system in- 
spired by conventional semiconductor doping. 

Semiconductors have narrow energy band gaps. A ma- 
terial may be intrinsically or extrinsically semiconduc- 
tive. Silicon, germanium, and gallium-arsenide are ex- 
amples of intrinsic semiconductors. Extrinsic semicon- 
ductors are doped with impurities. These impurities in- 
troduce allowed energy states within the band gaps of 
the host material. The introduced energy band gap, 
known as dopant-site bonding energy, is characteristi- 
cally small so that it takes little energy to ionize the 
dopant atoms and create free carriers in the conduction 
or valence bands. For example doping silicon with arsenic 
will create free electrons in the conduction band (n-type 
semiconductor). Doping with boron will create positive 
hole carriers in the valence band (p-type semiconductor). 
A typical level of impurity concentration would for exam- 
ple be 1 boron atom to 10 5 silicon atoms. 

Here we show how this property of semiconductors can 
be captured in a solid light system. Motivated by the 
semiconductor doping process, our analog process sim- 
ply involves changing the atomic transition of some of 
the two level atoms of the JC cavities. A tunable way 
of doing this is by Stark shift. We consider the ID case 
where every second site is detuned, but not far from res- 
onance as shown in Fig. [5] Although the 1:1 concentra- 
tion level is explicitly treated here, the physical proper- 
ties which we show here are general to other doping con- 
centrations. Embedding of a two-level system in waveg- 
uides [37j | and in a single resonator in coupled-resonator 
waveguides have been proposed for controlled scattering 
of photons [38]. Our model differs here in that all our 
cavities are embedded with two-level atoms, and we are 
controlling the detuning of multiple cavities. 

The unit cell, i.e. the minimum site grouping to form a 
periodic lattice, now constitutes two adjacent sites. Thus 
we need to extend the Bloch state of Eq. © to a two site 
basis, 

|</>,s)<g>|0,s+l) = \(f>,r}®\<l),r+l} exp{ik(d s -dr)], (19) 
and the onsite Hamiltonian becomes, 

2 

H r,r = ^{H 3 r c - uN r ) - K(a\a 2 + a x a\). (20) 

r=l 

The hopping Hamiltonian does not change. 

The Hamiltonian of Eq. (fTTj) will now be a (4n x An) 
block diagonal matrix, where the diagonal blocks are (4 x 
4) matrices. As before the blocks are the same function 
of n. Although an analytical solution is straightforward 
to obtain, it involves many terms and yields no insight, 




FIG. 5: (Color online) Every second atomic site is detuned 
from, but not far from, cavity resonance (Ao = 0, Ai = 
/3). The unit cell (dotted box) constitutes two adjacent sites. 
This introduces extra energy bands within the energy gap of 
the tuned system, in analogy to the semiconductor doping 
process. 

so we do not state it here. Instead we use it to plot the 
polariton-hole band structure and Mott lobes in Fig. [B] 
for tuned (A = A x = 0), doped (A = 0,A X = and 
detuned (Ao — Ai — (3) systems at different points of 
the phase diagram. 

Fig. E^a) shows the polariton-hole band structure for 
point I in Fig. Efd). The doped system has an extra 
polariton (hole) band. The extra band reduces the band 
gap from the tuned system. Here the lowest energy band 
is that of an extra polariton. This is analogous to the n- 
type semiconductor. Contrastingly Fig. EKb) shows the 
polariton-hole band structure for point II in Fig. [U(d), 
where the lowest energy band is that of a hole. This is 
analogous to the p-type semiconductor. 

Interestingly the solid light system also exhibits a fea- 
ture with no semiconductor analog. Fig. [6]Jc) shows the 
polariton-hole band structure for point III in Fig. [UJa). 
Here the doped system increases the band gap from both 
the tuned and detuned system, extending the MI phase 
of the doped system into the SF region of the tuned and 
detuned systems. This new effect speaks to the richness 
of the JCH systems and should open new avenues for 
investigation. 

The filled regions in Fig.[6]Jd)-(e) are the Mott lobes of 
the doped system. Because the detuning is not far from 
resonance, the Mott lobe states are: |n)A ® |ti)au|w ± 
l)A ®l n )Ai, where |n)A is the JC eigenstate |— , n) with 
Ao detuning. 

For a given n, the Mott lobe states will be, 

Ha ® Nai, |« + 1)a ® Had 

if n c (n, A ) < n c (n,Ai), (21) 
Ha ® \n)A 1 ,\n) Ao (g) |n + l>A 1 , 

if /i c (n,A ) > ft.(n, Ai), (22) 

In Fig. Efd)-(e) where A = and Ai = 0, the Mott 
lobe states are, 

|ti)a ® |n)Ai, |n + l)A ® Wah for n > 0, (23) 
|«)a ® Nai, |n)A ® \ n + 1)ai, for n = 0. (24) 
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FIG. 6: (Color online) (a) The solid (dashed) lines are the polariton (hole) band structure at point I. The left /middle/right 
panel is the band structure for the tuned/doped/detuned systen respectively. The impurity has introduced extra energy bands 
in the doped system reducing the energy band gap. The lowest energy band in the doped system is that of an extra polariton. 
This is analogous to a n-type semiconductor, (b) The solid (dashed) lines are the polariton (hole) band structure at point II. 
The lowest energy band in the doped system is that of a hole. This is analogous to a p-type semiconductor, (c) The solid 
(dashed) lines are the polariton (hole) band structure at point III. Interestingly the impurity increases the band gap: this has 
no semiconductor analogy, (d) The filled regions are the \n)A ® |«}ai (grey) and |n + 1)a ® I")ai (light grey) Mott lobes 
of the doped system, (e) The filled regions are the \n)A (8> |w)a! (grey) and \n — 1)a <S> |w)ai (dark grey) Mott lobes of the 
doped system. Overlayed in (d) and (e) in solid (dashed) lines is the phase transition border for the tuned (detuned) system. 
Note the different scales of (d) and (e). For notational convenience |ji)a ® |m-)Ai is written as |n)|m). 



In Fig.[Htd)-(e), overlayed in solid (dashed) lines is the 
Mott borders of the tuned (detuned) system. For small 
hopping the set of Mott lobes of the doped system (M 01 ) 
is the intersection of the set of Mott lobes of the tuned 
(M°) and detuned systems (M 1 ), M 01 = M°nM 1 , thus 
there are twice as many Mott lobes of smaller size in 
the doped system then in the tuned or detuned system. 
As hopping increases, the inter-cavity correlation causes 
divergence from this relation. Since E^ h (fi c ) ~ 0, \[i c — 
Mlmin gives the band gap. Hence the reduction in band 
gap which makes the Mott phase less stable to thermal 
fluctuations in the doped system is directly indicated by 
the smaller Mott lobe widths. 



the band structures we have shown the phase transition 
diagram, highlighting levels of MI stability. We have fur- 
ther demonstrated the similarities between the fields of 
solid state and solid light by showing a solid light analogy 
of the extrinsic semiconductor. Controlled site detuning 
can be seen as analogous to semiconductor doping, and 
we find that this process reduces the energy band gap 
forming smaller Mott lobes, but unlike the semiconduc- 
tor, it can also increase the energy band gap, extend- 
ing the Mott lobe tip. These results open the way for 
the treatment of other condensed matter systems in this 
framework. 
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